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We present an exact diagrammatic approach for the dimer-dimer scattering problem in two or 
three spatial dimensions, within the resonance approximation where these dimers are in a weakly 
bound resonant state. This approach is first applied to the calculation of the dimer-dimer scattering 
length as in three spatial dimensions, for dimers made of two fermions in a spin-singlet state, with 
corresponding scattering length ap and the already known result as ~ 0.60 hf is recovered exactly. 
Then we make use of our approach to obtain new results in two spatial dimensions for fermions 
^ ' as well as for bosons. Specifically, we calculate bound state energies for three bbb and four bbbb 

' resonantly interacting bosons in two dimensions. We consider also the case of resonant interaction 

between fermions and bosons and we obtain the exact bound state energies of two bosons plus one 
■ fermion bbf, two bosons plus two fermions bf^bfi, and three bosons plus one fermion bbbf. 

PACS numbers: 03.75.Ss, 05.30.Fk, 21.45.-l-v, 31.15.-p 
, Keywords: 
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I. INTRODUCTION 

Following the experimental realization of the Bose-Einstein condensation in ultracold bosonic gases, together with its 
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intensive study, the physics of ultracold Fermi gases has taken off recently with a strong development of experimental 
and theoretical investigations within the last few yearsi. In particular, much advantage has been taken of various 
■ Feshbach resonances which offer the possibility observing experimentally the so called BEC-BCS crossover. This has 
, been done in particular in ^Li and in '^''K. In the weak coupling limit of small negative scattering length, which is 
' realized far away on one side of the resonance, the corresponding weak attractive interaction between fermions leads 
10 to a BCS type condensate of Cooper pairs. On the other side of the resonance, where the scattering length is positive, 
weakly bound dimers, or molecules, consisting of two different fermions are formed. When one goes far enough of the 
resonance on this positive side, one obtains a weakly interacting gas of these dimers, which may in particular form a 
Bose-Einstein condensate, as it has been recently observed experimentalljiSiiii^. 

In the present paper, motivated by the problem raised by the physics of this dilute gas of composite bosons, we 
, will deal with the dimer-dimer elastic scattering and present an exact diagrammatic approach to its solution. This 
(— I ' will be done by staying in the so-called resonance approximation which is quite suited to the physical situation found 
Q [ with a Fesbach resonance. It this case the (positive) scattering length greatly exceeds the characteristic radius tq for 
O ' the attractive interaction between fermionic atoms. A problem of this kind was first investigated by Skorniakov and 
Ter-Martirosian& in the case of the 3-body fermionic problem. They showed that the scattering length of a fermion 
on a weakly bound dimer is determined by a single parameter, namely the two-body scattering length ap between 
fermions, and it is equal to l.lSap in the zero-range limit for the interatomic potential. A similar situation is found in 
the case of four fermions, where the dimer-dimer scattering length is fully determined by this same scattering length 

In a study of the crossover problem Haussmann^ calculated this scattering length of composite bosons as at the 
level of the Born approximation and found it equal to 2a p- This result was later on much improved by Fieri and 
StrinatiS., who took into account the repeated scattering of these composite bosons in the ladder approximation. 
This diagrammatic approach led them to a scattering length approximately equal to — O.lbap. However, this 
ladder approximation is not exact, because it misses an infinite number of other diagrams which in principle lead 
to a contribution of the same order of magnitude as those taken into account. Very recently this problem has been 
solved exactly by Petrov, Salomon, and Shlyapniko'\i2ii£ who found for the scattering length of these composite bosons 
aB = O-dap. This has been achieved by solving directly the Schrodinger equation for four fermions, using the well- 
known method of pseudopotentials. Here we will give an exact solution of this scattering problem of two weakly bound 
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T,(P) = ^:x^^xC 

FIG. 1: The graphic representation of the two particles vertex T2{P) (the four external propagators do not belong to T2{P)). 

dimers, using a diagrammatic approach in the resonance approximation, which can be seen as a bridge between the 
approach of Fieri and Strinat:^ and the exact result of Petrov, Salomon, and Shlyapniko\i2ii£. 

In order to show the strength and the versatility of our approach, we make use of it to obtain new results for 
various systems, in the two-dimensional (2D) case which is of interest not only for cold gases, but also for high 
superconductivity. Specifically we consider first a system of resonantly interacting bosons. We calculate exactly the 
three bosons hbb and four bosons hhhh bound state energies in this case. We also make use of our approach for the 
study of 2D bosons interacting resonantly with fermions. In this case we calculate exactly the bound state energies 
of two bosons plus one fermion hhf , two bosons plus two fermions bf^ bfi, and three bosons plus one fermion bbbf . In 
this respect the present paper is in the line of previous results obtained by of some of us. Indeed the possibility of 
two fermionsiiiiS // and two bosonsi^ bb pairing was predicted, as well as the creatioiiii of a composite fermion bf 
in resonantly interacting (a 3> tq) 2D Fermi-Bose mixtures. 



II. THREE PARTICLES SCATTERING 



As a preliminary exercise we will rederive the result of Skorniakov and Ter-Martirosian for the dimer-fermion 
scattering length as using the diagrammatic method^. Following Skorniakov and Ter-Martirosian, in the presence 
of the weakly bound resonance level —Ef, (with Eh > 0), we can limit ourselves to the zero-range interaction potential 
between fermions in the scattering of these two particles . The two-fermion vertex can be approximated by a simple 
one-pole structure, which reflects the presence of the s-wave resonance level in the spin-singlet state, and is essentially 
given by the scattering amplitude, namely: 

T2a/3;75(-P) = T2{P) X ((5„,^(5^j,5 ~ 5^, gS .{6a,^ S f) .1 + (5a.|(5^,|) = T2{P)x{a, f3)xh , ^) , (1) 

— 5171 — — ' xia,f3)^Sa,^^5p,i-Sa.iSp.T, (2) 
m'^'^ E — f'^/'lm + Eh 

where P = {P, E}, E is the total frequency and P is the total momentum of incoming particles, m is the fcrmionic 
mass, Eh — l/maj^. Indices a,/? and 7,5 denote the spin states of incoming and outgoing particles. The function 
x(q;, (3) stands for the spin singlet state. We will draw this vertex in the way, shown on Fig. ^ where the double line 
can be regarded as a propagating dimer. 

The simplest process that contributes to dimer-fermion interaction is the exchange of a fermion. We denote the 
corresponding vertex as A3 and it is described by the diagram on Fig. |5J Its analytical expression reads 

^3a,l3{Pl,P2;P) = -Sa,i3G{P - Pi - P2) , (3) 

where G{p) — 1/ {oj — p^/2m -I- ?0+) is the bare fermion Green's function. The minus sign in the right hand side of 
Eq.Q comes from the permutation of the two fermions. In order to obtain the full dimer-fermion scattering vertex T3 
we need to sum up all possible diagrams with indefinite number of A3 blocks. In the present case these diagrams have 
a ladder structure. It is obvious that the spin projection is conserved in every order in A3 and thus Tsa,i3 = ^a./S ?3. 
An equation for T3 will have the diagrammatic representation shown in Fig. It is obtained by writing that either 
the simplest exchange process occurs alone, or it is followed by any other process. In analytical form it reads 

nipi,P2;P) = -G{P - pi - P2) -Y,G{P - pi - q)G{q) T2{P - q) Ts{q,P2;P), (4) 

9 

where X] = qdVl / . We can integrate out the frequency Vl in Eq.Q by closing the integration contour in 
9 

the lower half-plane, since both T2{P — q) and Ts{q,P2\P) are analytical functions of SI in this region (this property 
for Ts{pi,p2\P) results from Eq.Q) itself). Hence only the "on the shell" value T-i{{q^, q^ / 2m] , p2] P) comes in the 
right- hand side of Eq.Q. Moreover, if we are interested in the low-energy s-wave dimer-fermion scattering length 
a3, we have to put P = {P, E] = {0, —Eh] and p2 = 0. Hence Eq.Q) reduces to an equation for the "on the shell" 
value of T3{pi,p2', P). Taking into account the standard relation between T-matrix and scattering amplitude (with 
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FIG. 2: The graphic representation of the simplest dimer-fermion scattering process A3 (the two external fermion propagators 
and the two external dimer propagators do not belong to A3). 




FIG. 3: The diagrammatic representation of the equation for the full dimer-fermion scattering vertex T3. 



reduced mass) and the fact that, from Eq.(^, T2 has an additional factor Sir /(rn^aF) compared to a standard boson 
propagator, we find that the fuU vertex T3 is connected with 03 by the following relation: 

^T3(0,0;{O,-£;fe}) = — 03. (5) 



This leads to introduce a new function 03 (k) defined by 
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a3(k) = — + V3W4+^j ({k, e/2m}, 0; {0, -E^}) . (6) 

and substituting it in Eq.Q, we obtain Skorniakov - Ter-Martirosian equation for the scattering amplitude: 



+ y/3k^/4 + mEb + mEb J {2itY {k^ + + k.q + mEb) 

Solving this equation one obtains the well known result^ for the dimer-fermion scattering length 03 = 03(0) = l.lSap 



III. DIMER - DIMER SCATTERING 



By now we can proceed to the problem of the dimer-dimer scattering. This problem was previously solved by 
Petrov et al^*iS via studying Schrodinger equation for a 4-fermions wave function. Our diagrammatic approach is 
conceptually close to Petrov's one. Its basic point is that it requires the introduction of a special vertex which describes 
an interaction of one dimer as a single object with the two fermions constituting the other dimer. 

Let us investigate all the possible types of diagrams that contribute to the dimer-dimer scattering vertex T^. In this 
process both dimers are temporarily "broken" in their fermionic components, which means that the fermions of one 
dimer exchange and/or interact with the fermions of the other dimer. The simplest process is an exchange of fermions 
by two dimers shown on Fig. ^Js. More complicated diagrams are composed by introducing intermediate interactions 
between exchanging fermions (see Fig. 01), c). As long as one of the fermions does not interact or exchange with the 
other ones, all these complications can be summed up in the T3 block (see Fig. 0Ji) which describes, as we have seen 
in the preceding section, the scattering of a fermion on a dimer. Furthermore we may exchange bachelor fermions 
participating in the T3 scattering. The resulting series has the diagrammatic structure shown on Fig. ^js. This series 
describes a "bare" interaction between dimers. The last obvious step is to compose ladder type diagrams from this 
" bare" interaction. A typical ladder diagram is shown on Fig. 0Jf. These general ladder diagrams describe all possible 
processes which contribute to the dimer-dimer scattering. 

The fact that the T4 vertex should be expressed in terms of T3 was first noticed by Weinberg in his work on 
multiparticle scattering problemsi&. Note that a calculation of the diagrams shown on Fig.0^, /requires information 
about an off-shell matrix T^, that is about a matrix with arbitrary relation between frequencies and momenta of 
incoming and outgoing particles. On the other hand, for the calculation of the dimer-fermion scattering length 03 in 
Eq.(I7|), only the simpler on-shell structure of T3 is required as we have seen in the preceding section. Luckily, as we 
will see now, we can exclude T3 from our considerations and express T4 only in terms of T2 . By doing this we reduce 
the number of integral equations required for the calculation of the dimer-dimer scattering length 04. 

Since, as we have just seen, it is impossible to construct a closed equation for the dimer-dimer scattering vertex 
T4, we wish to find an alternative way for taking into account in one equation all the diagrams contributing to 
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FIG. 4: The graphic representation of the dimer-dimer scattering processes contributing to T4,. 
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FIG. 5: Diagrammatic representation of the relation between the full dimer-dimer scattering matrix T4 and the vertex $. 

dimer-dimer scattering. Inspired by the work of Petrov et aliSiiS. and looking at the diagrams we have considered 
above, we are naturally lead to look for a special vertex that describes the interaction of two fermions, constituting 
the first dimer, with the second dimer taken as a single object. This vertex would be the sum of all diagrams with 
two fermions and one dimer as incoming lines. It would be natural to suppose that these diagrams should have 
the same set of outgoing - two fermionic and one dimer - lines. However in this case there will be a whole set of 
disconnected diagrams contributing to our sum that describe interaction of a dimer with only one fermion. As it 
was pointed out by Weinberg^^, one can construct a good integral equation of Lippmann-Schwinger type only for 
connected class of diagrams. Thus we are forced to pay our attention to the vertex ^apiQi,q2',P2, P) corresponding 
to the sum of all diagrams with one incoming dimer, two incoming fermionic lines and two outgoing dimer lines (see 
Fig.O. This is also quite natural from our view point since, in our scattering problem we are interested in a final state 
with two outcoming dimers. Indeed once this vertex ^afsili, 12', P2, P) is known, it is straightforward to calculate the 
dimer-dimer scattering vertex T4(pi,p2; P) which is given by: 

T4(pi,p2;P) = ^ x{a,m{P+Pi-k)G{k)^^p{P + pi-k,k;p2,P). (8) 

fe; Q,/3 

The corresponding diagrammatic representation is given in Fig. [5| One can readily verify that, in any order of 
interaction, $ contains only connected diagrams. 

The spin part of the vertex $a,/3 has the simple form ^a,p{<lii<l2] PtP2) — xic^i l3)^iQi7<l2', P,P2)- The diagrammatic 
representation of the equation for 4> is given in Fig. |S1 One can assign some " physical meaning" to the processes 
described by these diagrams. The diagram of Fig. represents the simplest exchange process in a dimer-dimer 
interaction. The diagram of Fig. |3) accounts for a more complicated nature of a "bare" dimer-dimer interaction. 
Finally the diagram of Fig. allows for a multiple dimer-dimer scattering via a "bare" interaction (it generates 
ladder- type diagrams analogous to those of Fig. The last term in Fig. El means that we should add another set of 
three diagrams analogous to those of Fig. b, c but with the two incoming fermions (qi and 92) exchanged. The 
diagrammatic representation translates into the following analytical equation for the vertex $: 

$(gi, 92;P2, P) = -G{P -qi+ P2)G{P -q2-P2)-J2 G{k)G{2P - <zi - 92 - k)T2{2P - qi - fc)$(gi, k;p2,P) 

k 

-\Y.G(Q- 9i)G(2P - Q - q2)T2{2P - Q)T2{Q)G{k)G{Q - fc)$(fc, Q - k; p2, P) + (gi ^ 92). (9) 

Q,k 

Finally let us also indicate that it is possible to rederive the same set of equations, purely algebraically, by taking 
a complementary point of view. Instead of focusing, as we have done, on the free fermions lines as soon as a dimer 
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FIG. 6: The diagrammatic representation of the integral equation for the function $ introduced by dimer-dimer scattering . 

is "broken", we can rather keep track of the fermions which make up a dimer. This leads again automatically to 
introduce the vertex ^{qi,q2;p2, P)- Then Eq.® is recovered when one keeps in mind that, after breaking dimers, 
one may have propagation of a single dimer and two free fermions before another break (this corresponds to the second 
term in the right hand side of Eq.Q). Alternatively one may also have the propagation of two dimers, which leads 
to the third term in Eq.®. 

Coming back now more specifically to our problem, wc can put P2 ~ and P ~ {0, —Eb\ since we are looking for an 
s-wave scattering length. At this point we have a single closed equation for the vertex $ in momentum representation, 
which we believe is analogous to Petrov ei al equation in coordinate representation. To make this analogy more 
prominent we have to exclude frequencies from the equation by integrating them out. However this exclusion requires 
some more technical mathematics and we leave it out for Appendix A. 

The dimer-dimer scattering length is directly related to the full symmetrized vertex T4{pi,p2', P). Just as in the 
preceding section, taking also statistics into account, we have: 



If one skips the second term in Eq. (Q, i.e. one omits diagram Fig|3), one will arrive at the ladder approximation of 
Fieri and Strinati*. The exact equation ^ corresponds to the summation of all diagrams. We have calculated the 
scattering length in the ladder approximation and the scattering length derived from the exact equation and obtained 
O.lSap and 0.60ap respectively. Some details on our actual procedure are given in the next section. Thus our results 
in the ladder approximation are in agreement with the results® of Fieri and Strinati and, in the general form, with 
the results of Fetrov et. a/S^iS. Note also that our approach allows one to find the dimer-dimer scattering length in 
the 2D case (this problem was previously solved by Petrov et. aJi^). 

Finally we would like to mention that our results allow one to find a fermionic Green's function, chemical potential 
and sound velocity as a function of of in the case of dilute superfluid bose gas of dimers at low temperatures. The 
problem of dilute superfluid bose gas of di-fermionic molecules was solved by FopoAii^, and later deeply investigated 
by Keldysh and Kozlov^^. Those authors managed to reduce the gas problem to a dimer-dimer scattering problem in 
vacuum, but were unable to express the dimer-dimer scattering amplitude in a single two-fermion parameter. A direct 
combination of our results with those ones of Popov, Keldysh and Kozlov allows one to get all the thermodynamical 
values of a dilute superfluid resonance gas of composite bosons. Another interesting subject for the application of our 
results will be a high-temperature expansion for the thermodynamical potential and sound velocity in the temperature 
region T ^ ^ Eb, where the composite bosons begin to appear. 



Let us give now some details on the way in which we have solved effectively the above equations. Actually we 
have dealed with two problems, the scattering length calculation discussed above and the bound states problem to be 
discussed below. Our two problems are quite closely related since, for the scattering length problem, we look for the 
scattering amplitude at zero outgoing wavevectors and energy for two dimers, while for the bound states we look for 
divergences of this same scattering amplitude at negative energy. As already indicated, in both cases the situation is 




(10) 



IV. PRACTICAL IMPLEMENTATION 
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somewhat simplified with respect to the variables we have to consider, due to the specific problem we handle. First 
with respect to P = {P, E}, we have P = since we work naturally in the rest frame of the four particles. Moreover, 
with respect to the total energy, E = —eEi, is negative. Specifically e = 1 when we look for the scattering length. 
Or when we consider bound states e gives the energy of the bound states we are looking for. Next, with respect to 
parameter p2 = {p2jP2} which characterizes the outgoing dimers, we will have naturally p2 = as we have said 
since we consider zero outgoing wavevectors. Since we will evaluate p2 on the shell, we have merely p2 — 0, and this 
parameter drops out. Hence in the following we do not write anymore explicitely the value of parameter P. Both for 
the scattering length problem and the bound states problem, we have followed two main routes. 

In our first route, we have written a specific integral equation for T4{pi,p2), which is then solved numerically. The 
details of our derivation for this integral equation are given in Appendix B. The kernel for this equation is itself 
obtained from a vertex F. The defining integral equation Ea. (jB2|) for this vertex has been inverted numerically, 
by calculating the inverse matrix, to obtain the vertex F(gi,g2;P2)- We have used~, LU factorization and Gauss 
quadrature. The result has then been substituted in Ea. (|Bll) which gives the kernel A4(pi,p2) coming in the integral 
equation Ea. ljB3|l . The solution of this last equation is naturally also handled numerically, for example by finding the 
eigenvalues of the kernel for the bound states problem. 

In our second route we have kept both functions T4 and $. In the following we do not write anymore the parameter 
P2 which takes always the trivial value p2 — 0, as explained above. Hence we are left with T4{pi) which, because of 
rotational invariance, depends only on the energy pi and the modulus |pi| of the momentum. For brevity we denote 
this quantity t4(|pi|,pi). On the other hand it is shown in Appendix A that, in order to evaluate the second term 
in the right-hand side of Eq.©, we need only the evaluation of <i>(qi, (72) on the shell, which we denote as 0(qi, q2). 
It depends only on the three variables |qi|, |q2| and the angle between these two vectors. Hence it is enough to 
write Eq.Q only for qi and q2 taking on the shell values. From Eq.@ this leads for (/'(qi, q2) to the following more 
convenient equation: 

,^ . 1 f d^k 2m^2(2|g| + [2k2 + 2qH(k + qi)^]/4m)0(qi,k) 

n^^^'i^) {\E\+<il/m)i\E\ + qym)^ J (27r)^ 4to|S| + k2 + q2 + q2 + (k + qi + q2)2 ^ > 

1^ (2m)2^2(|^|+Q + QV4m)^2(|^|-Q + QV4m)t4(|Q|,Q) . . 

2 ^ {2m{\E\ _ Q) + q2 + (Q + qi)2)(2m(|i?| + Q) + + {Q + q^)^) + ^"^^ 

where the dimer propagator t2{x) depends on the space dimension D. For Z? = 3 it is given from Eq.l^ by t2{x) = 
—4:TT/[m^^'^{^/x— ^/Eb)], while for D = 2 according to Ea. (|^ we have t2ix) = —4'K/[mln{x/Eb)]. In the third term 
the angular integration can be performed analytically, and one is left with double integrals for the last two terms, for 
the 3D as well as for the 2D case. It is actually quite convenient, in the last term, to deform the Q contour from 
] — cxD, oo[ to ] — ioo, ioo[ by rotating it by n/2. No singularity is met in this deformation, and one is left to deal only 
with real quantities. 

The above equation has to be supplemented by a corresponding equation for t4(|q|, 9) obtained from the definition 
Eq.®. The important point is that the additional integrations can be performed analytically, owing to the various 
invariances under rotations found in the resulting terms. We just give here as an intermediate step the structure of 
the resulting equation: 

/•oo /"Oo /'27r 

h{k,iz)^ S{k,z)+ dpi dp2 rfa/(fc,z,pi,p2,a)i2(2|S| + [3pf + 3pl + 2pi.p2)2]/4m) ?!)(pi,p2) (12) 
Jo Jo Jo 



/>00 />CXD 

dK dZ J{k,z,K,Z)\t2{\E\+iZ + K^/Am)]^ t4{K,iZ) 
Jo Jo 



where a is the angle between pi and p2. Here S{k,z), I{k, z,pi,p2,a) and J{k, z, K, Z) are analytically known 
functions of the variables (except that J requires to perform numerically a simple integration to be obtained, see 
below). In this equation and in particular in its last term, we have already gone to the purely imaginary frequency 
variable for 14. The resulting t4{x,iz) turns out to be real and even with respect to z. 

To be fully specific let us now give the actual self-contained integral equations which we have solved. We restrict 
ourselves to the 3D case and to the bfbf case (implying a = 1), corresponding to the dimer scattering problem 
treated^iili by Petrov et al. The only generalization is that we keep E — — e|i?h|, instead of setting e = 1 as we should 
if we considered only the scattering length problem. For clarity we write the resulting equations with dimensionless 
quantities, where 1/a has been taken as unit wavevector, and \Eb\ = 1/ma? as energy unit. For simplicity we keep 
basically the same notations for the various variables. We just indicate by a bar over the function name that they 
are expressed in reduced units, with reduced variables (actually we write i4{k,z) instead of ^4(^,12), and there is a 
change of sign between 0(qi,q2) and (/>(pi,P2))- Equations for other cases and dimensions are completely similar 
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with only few changes in coefficients, signs (for the particle statistics), for the expression of t2{x) and for the explicit 
functions coming from analytical angular integrations. 
We obtain: 



, , 1 , 1 r,.2„. r . ^^^ 0(pi,k)f2(2e+(3pf + 3fc^ + 2fcpiCosg)/4) , ^ 



dz dkU{k,z)\t2{e + k^lA + iz)\^ I{Bi,B2,a) (13) 

T^PlP2 Jo Jo 

with A± = 2e + pI + P2 + k'^ +P1P2 cos a + kpi cos 9 + kp2 cos(q! ± 9), and a is the angle between pi and P2, while 9 
is the polar angle of k with pi. We have simply set now ^2(2;) = [1 ^ V^]^'^- Here we have also defined the function: 

/ Si,B2,a =ine^^ln— ^ 14 

2/B BlB*+cosa-^/E 

E ^ B'l + Bf + 2BiB^ cos a - sin^ a (15) 

B{p,k,z)^^[e + p^ + ^-iz] (16) 

B, = B{p,,k,z) (17) 

The corresponding equation for ti{k, z) is; 



- , 1 , 1 + cos 7 + 2^cos7COs((p — 7/2) 



47rfcz 1 + cos 7 + 2^cos7cos(iy9 + 7/2) 

-77^ / Pidpi / P2dp2 / sinada ^(pi,p2)t2(2e + (3pi + 3p2 + 2piP2COsa)/4)/(Bi,B2,a) 
1^ Jo Jo Jo 



1 



00 



KdK / dZf4(if,^)|i2(e + if74 + iZ)|V(fc,z,X,Z) (18) 
"'0 



27r3fc 

with (/9 = arctan(fc/2) and 7 = arctan[4z/(4 + fc^)] and we have defined the function: 

7(1- 7\ 1 g(^,fc,^,^) , C{x,K,k,z) 

J{k,z,K,Z)= dx p-^lnp-^ K- 7\ ^^7^( K- h \ 

Jo e + .x2 + -i-ti^ C{x,-k,K,Z) C{x,-K,k,z) 

C{x,k,K.,Z) = [e+{x + ^f + ^f + Z^ (20) 

It is seen on these integral equations for our two unknown functions ti{x,z) and 0(pi,P2) that they require only 
at most a triple integrals to be performed numerically. In this sense they are not numerically more complicated than 
the work involved in solving directly for the corresponding Schrodinger equation, as it has been don^ii^ by Petrov 
et al. Indeed these integrals require only a few appropriate change of variables to take care of singular behaviours 
occuring on some boundaries. Otherwise they have been performed with unsophisticated integration routine. 

In the case of the scattering length a mere iteration algorithm has been found to lead rapidly to the solution 
(provided an appropriate exact algebraic manipulation is made to make the iteration convergent). In this way we 
have been able to handle 45 x 45 x 45 matrices (for the three variables entering ^(pi, P2)). This size is large enough to 
allow improved precision by extrapolation to infinite size, although we have not done it in the present case, but rather 
for the ground state of the bhhh complex discussed below. This leads to the result as = 0.60 aj? in full agreement 
with Petrov^' et al, within a quite reasonable computing time on (nowadays) unsophisticated computer. We have 
not tried to improve on the accuracy of the result, since there is no basic interest. In the case of the bound states, to 
be described below, we have proceeded to a straight diagonalization of the matrix equivalent to the right hand sides 
of Ea. H13|) and Ea. H18|l with the Lapack library algorithm. In the 2D case, it is worth noticing that, because of the 
logarithmic dependence of ^2(2;) on x, it is quite an improvement to make the change of variables K = e^^'^K' and 
Z = eZ' , and so on, since the more appropriate variable turns out to be Ine rather than e itself. 



V. NEW RESULTS IN A 2D CASE 



We will now apply the diagrammatic approach developed in the previous sections (see also Appendix 0) to get new 
results for the systems of resonantly interacting particles in a 2D case. 
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As it was first shown by Danilo\i2i (see also a paper by Minlos and Fadee\i2S) in the 3D case, the problem of three 
resonantly interacting bosons could not be solved in the resonance approximation. This statement stems from the fact 
that in the case of identical bosons the homogeneous part of Skorniakov-Ter-Martirosian equation (jjj) has a non-zero 
solution at any energies. The physical meaning of this mathematical feature was elucidated by Efimov, who showed 
that a two-particle interaction leads to the appearance of an attractive 1 /r^ interaction in a three-body system. Since 
in the attractive potential a particle can fall into the center, the short range physics is important and one can 
not replace the exact pair interaction by its resonance approximation. 

On the contrary in the case of the 2D problem the phenomena of the particle fall into the center is absent and 
one can utilize the resonant approximationi22i24. Therefore it is possible to describe three- and four-particle processes 
in terms of the two-particle binding energy Et, = 1/ ma?' only (below, for simplicity we will assume that all particles 
under consideration have the same mass m) . We will leave aside the problem of composite particles scattering and 
will concentrate on the problem of binding energies of complexes of three and four particles. 

As well as in the case of the 3D problem, the cornerstone in the diagrammatic technique is the two-particle resonance 
scattering vertex T2 (see Fignj. For two resonantly interacting particles with total mass 2m it reads in 2D: 

47r a 

' " ^~^\i-i{{V^/Am^ E}l\EB\y ^ 

where we introduce a factor a = {1,2} in order to take into account whether two particles are indistinguishable or 
not. That is a = 2 for the case of a resonance interaction between identical bosons, while a = 1 for the case of a 
resonance interaction between fermion and boson, or for the case of two distinguishable fermions. 



A. Three particles in 2D 

We start with a system of three resonantly interacting identical bosons - bbb - in 2D. An equation for the dimer- 
boson scattering vertex T3 which describes interaction of three bosons has the same diagrammatic form as the one 
shown on the Figl^l however there are small changes in the rules for its analytical evaluation. The resulting equation 
reads: 

T3(pi,P2; P) = G{P -Pi-P2)+Y^ G{P - - q)G{q) T2{P - q) T^{q,P2; P), (22) 

q 

where we have now = * J (P qdil / {2tt)^ , P — {0,E}, and one should put a — 2 for the two-particle vertex T2 

in Eq. H21() . The opposite signs in Eq. for fermions and Eq. (|22() for bosons are due to the permutational properties 
of the involved particles : an exchange of fermions (see Fig|21) results in a minus sign, while an analogous exchange of 
bosons brings no extra minus. Finally, as we mentionned above, we note that three-particle s-wave (s-wave channel 
of a boson-dimer scattering) binding energies E3 correspond to the poles of T3(0, 0; {0, — 1£^3|}) and, consequently, at 
energies E = E3 the homogeneous part of Eq. (|22() has a non-zero solution. Solving Eq. (|22|l we find that a complex 
of three identical bosons has two s-wave bound states E^ = — 16.5i?f, and £'3 — — 1.27£'b in accordance with the 
previous results of Bruch and Tjon2i24. 

Let us now consider a complex - fbb - consisting of one fermion and two bosons. As noted above we take bosons 
and fermions with equal masses rrib = rrif = m. We assume that a fermion-boson interaction Ufb, characterized by 
the length r/f,, yields a resonant two-body bound state with an energy E = —Eb. In the same time a boson-boson 
interaction Ubb, characterized by the interaction length rbb, does not yield a resonance. Hence if we are interested in the 
low-energy physics the only relevant interaction is Ufb and we can ignore the boson-boson interaction Ubb, the latter 
would give small corrections of the order |£^B|mr^j <C 1 at low energies. In order to determine three-particle bound 
states one has to find poles in the dimer-boson scattering vertex T3. Since we neglect the boson-boson interaction 
Ubb the vertex T3 is described by the same diagrammatic equation of Fig. |31as for the problems of three bosons. The 
analytical form of this equation also coincides with Eq. H22() with the minor difference that the resonance scattering 
vertex T2 now corresponds to the interaction between a boson and a fermion, and therefore we should put a = 1 in 
Eq. (|21(l for T2. Solving the equation for T3 we find that fbb complex has only one s-wave bound state with the energy 
E3 = —2.39 Eb- Note that a complex - bff - consisting of a boson and two spinless identical fermions with resonance 
interaction Ufb does not have any three-particle bound states. 



B. Four particles in 2D 



After solving the above three-particle problems we may proceed to the complexes consisting of four particles. At 
first we will consider four identical resonantly interacting bosons bbbb^^. Any two bosons would form a stable dimer 
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with binding energy E = —E^. We are going to find a four-particle binding energy as an energy of an s-wave bound 
state of two dimers. Generally speaking a bound state could emerge in channels with larger orbital moments, however 
this question will be a subject of further investigations. Just as in the preceding subsection, in order to find a binding 
energy we should examine the analytical structure of the dimer-dimer scattering vertex T4 and find its poles. The 
set of equations for T4 has the same diagrammatic structure as those shown on Fig. O and Fig. The analytical 
expression for the first equation reads: 

and the equation for the vertex $ is: 
$(gi,g2;P2,P) =G(P-gi+p2)G(P 

" Q,k 

(24) 

where T2 should be taken from Ea. H21() and one should put a = 2 for the case of identical resonantly interacting 
bosons. When we look for the poles of T4 as a function of the variable E, with P — {0, i?}, we have naturally to 
consider only the homogeneous part of this equation. We have found 2 bound states for the &66fe complex. The values 
of the total binding energy |i?4| — 2\E\ are given in Table 1 below. Certainly for the validity of our approximation we 
should have |i?4| <C l/mrg. For the case of four bosons hhbh it means that 197 ^ l/mrg and hence a/rp ^ -\/T97. 
This case can still be considered as quite realistic for the Feshbach resonance situation. 

The case of a four-particle complex - bf-^bfi - consisting of resonantly interacting bosons and fermions is still 
described by the same equations H23I24|) but with parameter a — 1. In this case we found 2 bound states and they 
are also listed in Table 1. 

In order to obtain bound states of the fbbb complex one has to find energies P — {0, E} corresponding to nontrivial 
solutions of the following homogeneous equation 

$(qi, 92; P2, P) = Y. G{k)G{2P -qi-q2-k) T2{2P - qi - k) k-p2,P) + {qi ^ 92). (25) 

k 

This equation corresponds to the diagram of Fig. Eb- We have found a single bound state for this fbbb complex. 
Finally wc summarize the results concerning binding energies of three and four resonantly interacting particles in 2D 
in Table 1. 



Table 1. Bound states of resonantly interacting particles in 2D 



System 


Relative^^ 
interaction 


Number of 
bound states 


Energy (in \Eb\)^^ 




bbb 


Ubb 


2 


1.27 , 16.5 


2 


fbb 


Ufb 


1 


2.39 


1 


fbbb 


Ufb 


1 


4.1 


1 


bfMl 


Ufb 


2 


2.8, 10.6 


1 


bbbb 


Ubb 


2 


22. , 197. 


2 



Interaction that yields resonance scattering. All other interactions are negligible. 



m = mi, = mj. 
■^The indistinguishability parameter in Eq. 1211 . 

For the bbbb complex we find the beginning of a contimmm of states at \E4\/Eb — 16.5, as it should be since 
this is, within our numerical precision, the binding energy of bbb. Similarly we find the beginning of a continuum 
at lEil/Eb = 2.4 for the fbbb and the bf^jbfi complex, in agreement with the binding energy of fbb. We display 
our corresponding results in Fig. [7| and Fig.|Sl In all our calculations we find numerically, as a function of |i?4|, the 
eigenvalues A corresponding to the matrix on the right-hand side of our equations, for example Eq. H25|l . When one 
of these eigenvalues is equal to 1, this means that the corresponding E4 is the energy of an eigenstate of our complex. 
In Fig. [7| we display the first highest eigenvalues for \E4\ = 2.4, both for the bf^bf^ case and the bbbf case. One sees 
clearly that a fair number of eigenvalues are essentially equal to 1 . One could tune them exactly to 1 by changing very 



- V G(P +pi - fc)G(fc)$(P +pi-k, k;p2,P), (23) 



" k 



-q2~P2)+Y, G(fc)G(2P -qi-q2- fc)T2(2P - - fc)$((?i, fc;p2, P) 

k 

- q2)T2{2P - Q)T2{Q)G{k)G{Q - k)<^{k, Q - k; p2, P) + (gi ^ 92). 
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3,5 
3 
2,5 

e< 2 

1,5 
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0,5 
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FIG. 7: Eigenvalues A found for |i54| = 2.4, both for the bbhf case and the bf^bfy case. When one of the eigenvalues is equal 
to 1, Ei is the energy of a possible eigenstate of the complex. The number n appearing on the x axis is just here to number 
the first few eigenvalues which are displayed by decreasing order. 
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FIG. 8: Eigenvalues A found for |i54| = 16.5 and |i54| = 22., for the bbbb case. When one of the eigenvalues is equal to 1, Ei is 
the energy of a possible eigenstate of the complex. The number n appearing on the x axis is just here to number the first few 
eigenvalues which are displayed by decreasing order. 

slightly |i?4|. Hence this corresponds to the beginning of the continuum. By contrast one sees also clearly two isolated 
eigenvalues larger than 1, for the bf-^bf ^ case, and one eigenvalue larger than 1 for the bbbf case. One can bring them 
to A = 1 by increasing |i?4|, and therefore they correspond to the bound states that we have found. Similarly we 
display in Fig. |51the eigenvalues for the bbbb case, for the value \E4\ = 16.5 corresponding essentially to the threshold 
for the continuum. Here again one sees many eigenvalues quite close to 1. On the same figure we also show the results 
of the same calculations for \Ei \ = 22. in order to display the way in which this whole spectrum evolves with |i?4|. In 
particular one sees clearly the two isolated eigenvalues, corresponding to the two bound states found in this case. In 
particular since one of them is equal to 1, this means that the binding energy of one of the bound states is equal to 
22 Eb, within our numerical precision. 

Note finally that all our calculations correspond to the case of particles with equal masses mf = nib = although 
they can be quite easily generalized to the case of different masses. 
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VI. CONCLUSIONS AND DISCUSSION 



For the problem of resonantly interacting fermions in 3D we have developed an exact diagrammatic approach that 
allows to find the dimer-dimer scattering length ag — 0.600^^ in exact agreement with known results. This exact 
diagrammatic solution of the dimer-dimer scattering length problem in 3D opens new horizons for the extension of 
the self-consistent mean-field schemes of Leggett and Nozieres-Schmitt-Rink to the inclusion of the quite essential 
three and four-particle physics in the two-particle variational wave-functions of the BCS-type. This in turn will help 
us to get diagrammatically exact results for Tc, pseudogap and sound velocity in the dilute BEC- limit and to develop 
a more sophisticated interpolation scheme for these quantities toward the unitarity limit. The work on this very 
exciting project is now in progress. 

We have applied the developed approach to get new results in the 2D case. Namely, we have calculated exactly the 
binding energies of the following complexes: three bosons bbb, two bosons plus one fcrmion bbf , three bosons plus one 
fermion bbbf, two bosons plus two fermions bf^bfi, and four bosons bbbb. 

Our investigations enrich the phase-diagram for ultracold Fermi-Bose gases with resonant interaction. They serve as 
an important step for future calculations of the thermodynamical properties and the spectrum of collective excitations 
in different temperature and density regimes, in particular in the superfluid domain. Note that in purely bosonic 
models in 2D or in the Fermi-Bose mixtures in the case of prevailing density of bosons ng > np a. creation of larger 
complexes consisting of 5, 6 and so on particles is also possible. In fact here we are dealing with the macroscopic 
phase separation (with the creation of large droplets). The radius of this droplet Rn ioi N bosons in 2D is estimated 
iripSk on the basis of a variational approach. Note that already for = 5 the exact calculation of the bound state 
requires huge computational capability, but it would be interesting to see precisely how this would appear with our 
approach. 



APPENDIX A: DIMER-DIMER SCATTERING EQUATION. FREQUENCY INTEGRATION 

In this Appendix we will show how one can integrate explicitely over the frequency dependence in the dimer-dimer 
scattering equation ^ (we consider only this case, the other ones considered in Section IV would require trivial 
modifications). To simplify further computations we slightly change the notation and introduce a chemical potential 
fi = —Eb/2 and the single fermion energy = p^/2to — h = p^/2m -f -Efc/2, with the modified fermion Green's 
function Q{p) = l/{uj ~ S^p). In the expression Eq.QJ for 72(Q) we have similarly to replace E hy E — Ei,. 

The integral equation ^ reads more explicitely (with k = {k, lj} and Q = {Q, il}): 

oo 

<?(9l,92) = -5(-gi)e(-g2)-* J ^ J -^gik)gi~q,-q2-k)T2{-qi-k)<l>{quk) + 

— OO 

+ IJ ^-0^^Q{Q~qi)Q{^Q-q2)r2{-Q)r2{Q)g{k)g{Q~k)<P{k,Q~-k) + {q,^q2). (Al) 

From this equation <?(gi, (72) — ^(<Z2, 91), as it is obvious physically. Note also that the third term is already explicitely 
symmetrical in qi <~> (72- 

First we note that, from Ea. (|All) itself, <f (51,(72) is analytical with respect to the frequency variables loi and W2 
of the four- vectors qi and (72 in the lower half-planes ^vaui < and 3m W2 < 0. This can be seen by assuming this 
property self-consistently in the right-hand side, and checking that the three terms are then indeed analytical, or 
equivalently one can proceed to a perturbative expansion. Then, if we make the "on the shell" calculation of 'l>{qi, 92) 
from Ea. (|Aip . that is for ui — and llJ2 — ^q2: we see that, for second term in the right-hand side, the only 
singularity in the lower complex plane Jmuj < is the pole of g{k) at lu ~ ^k- Hence the integration contour can be 
closed in the lower half-plane, leading to: 

00 

^ f ^gik)g{-q,-q2-k)U-q,-k)Hq,,k)^~ j>(qi,k). (A2) 

— 00 

Here we denote <?(qi,q2) = '^({qi, CqJ, {q2, CqsD- 

The frequency integration of the third term in Ea. (jAl(l over the frequencies fl and cj is more difficult because 
singularities are not essentially located in one half of the complex plane, as it was the case for the second term. For 
example <l>(k,Q — k) has singularities in both half planes, with respect to w, and similarly for 72(— Q)72((5) with 
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respect to il. We solve this problem by splitting the involved functions as the sum of two parts, one analytical in the 
upper complex plane, and the other one in the lower complex plane. 

First we write F{Q, Q, Qi, = G{Q — qi)G{—Q — <l2)T^2{—Q)T^2{Q) + (^i ^ 92) (we take into account that we want 
to calculate (?2)"on the shell") as: 

F{^, Q, Qi, Qa) = U+{n, Q, qi, qa) + U^{^1, Q, q^, q^) (A3) 

where 11+ and U- are respectively analytical in the upper and lower complex planes of 17. This is done by making 
use of the Cauchy formula /(fi) = (l/2i7r) dz f{z)/{z — Vl) for a contour C which encircles the real axis (on which 
F has no singularity) and is infinitesimally near of it. This gives: 

Mf^,Q,q,,q,).-L r,.£lflQlSlZ^ (A4) 

with e = 0+. Making use oi F{-^) = F{n), we find (7_(r2, Q, qi, qa) = ;7+(-0, Q, q^, qj). 

On the other hand the last part of the third term T4{Q') = J d^k' g{k')g{Q' - k') <P{k' , Q' - k') = J d^k' g{Q' /2 + 
k')g{Q'/2 - k')${Q'/2 + k',Q'/2 - k') satisfies T_i{-Q') = Ti{Q'). This can be seen by substituting Eq.(|JlJ for 
<1>{Q' /2 + fc', Q' /2 — k') in this last expression for Ti{Q'). For the first term contribution, the result is trivial. For the 
second term, one has to make the shift fc fc — Q'/2, and then fc <-> fc'. In the third term one has to make the shift 
k' ^ k' + Q/2 and then k' —k'. Then, when we make the change Q — Q in the third term of Ea. llAl|) and use 
Fii—Q) = Ti{Q)^ we see that the C/_ contribution is exactly identical to the U+ contribution and we are left with a 
single contribution from L/_ to evaluate. 

In order to perform the uj integration in f^iQ) = Jd'^kg{Q/2 + k)g{Q/2 - k)${Q/2 + k, Q/2 - fc), we split: 

'P{Q/2 + k, Q/2 -k) = 'P+iQ/2 + fc, Q/2 - k) + 'P-{Q/2 + k, Q/2 - k) (A5) 

into the sum of two functions, with ^+ analytical in the upper complex plane with respect to u>, and analytical 
in the lower complex plane. That this can be done is immediately seen from Ea. (|Al|l itself. For the first term we just 
have to write the product of Green's functions as g{k - Q/2)g{-k - Q/2) = ~{g{k ~ Q/2) + g{-k - Q/2))/{^l + 
Ck+Q/2+4(k-Q/2 — *e): which has explicitely the required property. In the third term we can handle the product of the 
first two Green's functions in the same way. Finally, in the second term, after performing the uo integration as indicated 
above (but without taking the " on the shell" values for the frequencies) , one sees that the result for the term written 
explicitely above in Ea. ljAl|) is analytical in the lower complex plane with respect to uj. The corresponding term 
obtained by (qi (72) is analytic in the upper complex plane. In each case one checks that the functions analytical 
in the upper and lower complex plane are related by fc — fc, so that <1>- {Q/2 + k,Q/2 — k) = <!>+{Q/2 — fc, Q/2 + k). 
Hence by the change of variable k ^ —k, the contributions of and are equal. 

Then we have arrived, for the calculation of 14 (Q), to a situation which is similar to the one we met for three 
particles. Since <P+{Q/2 — k,Q/2 + k) and g{Q/2 — k) are analytical in the lower complex plane, we can close the 
integration contour at infinity in this lower half plane and the only contribution comes from the pole of g{Q /2 + k). 
This leads to: 

J (27r)-^ il - 4k+Q/2 - 4k-Q/2 + «e 

where k, Q) is 'l>+{Q/2 — k,Q/2 + k) evaluated for uj = Ck+Q/2 ^ fl/2. An important property, which can be 
checked on each term contributing to <P+{Q/2 — fc, (5/2 + k) is that J-"(r2,k, Q) is analytical in the lower complex 
plane with respect to Vl. Hence the integration of C/_(ri, Q, q^, q2)T4((5) over can also be performed by closing 
the contour in the lower half plane, since the only singularity in this half plane is the pole due to the denominator 
in Ea. (|A6p . The contribution of this pole leads to the evaluation of JF(ri, k, Q) for Q, = '^k+Q/2 + Ck-Q/2- Taken 
with the above definition of T this means that we have calculated <l>+{Q/2 — /c, Q/2 + k) for Vl/2 — ui ~ Ck-Q/2 and 
il/2 + u! = Ck+Q/2: which is just an evaluation "on the shell". Because of the simple relation between and the 
result can be expressed in terms of <?(k + Q/2, — k + Q/2) itself. 

Gathering all the above results we end up with the following complete equation for <?(qi, q2): 

^ - I ^2 (-Cqi - Ck, qi + k) 

4Cqi4q2 J i^T^r ?qi + ?q2 + 4k + Cqi+q2+k 

d^Q d^k 



/d'^O d"^k 
^ (^C^(^k+Q/2 + ^k-Q/2, Q, qi, q2) <Z'(k + Q/2, -k + Q/2) + (q^ ^ q^). (A7) 
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FIG. 9: The diagrammatic representation of the equation for the full dimer-dimer scattering vertex T4{pi,p2; P). 
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FIG. 10: The diagrammatic representation of the sum of all irreducible diagrams A4{pi,p2; P). 



In this equation we have modified the integration contour in the definition of U- to have it running on the imaginary 
axis rather than on the real axis, and we have used the symmetry property of _F(z, Q, qj^, qj) with respect to z, 
together with symmetry properties of ??(qi, q2), to rewrite the result in terms of the real function: 



C/(17,Q,qi,q2) 



n 



dy 



■Ffa,Q,qi,q2) 



(A8) 



which shows that <?(qi, q2) itself is real. 

We have made practical numerical use of Eq. l|A7p to find for example the ground state energy. Although this 
turned out to be quite feasable, this equation appears finally less convenient than what we have described in section 
IIVI This was expected since the solution implies quadruple integrals, instead of the triple integrals we had only to 
deal with in section Hvl 



APPENDIX B: MODIFIED DIMER-DIMER SCATTERING EQUATION 

This appendix is devoted to an alternative description of the dimer-dimer scattering process. The purpose is to 
obtain a direct integral equation for T4^{pi,p2; P), in a way convenient for numerical calculations. Below we derive 
such a set of equations, that were used for practical computations as indicated in section HVI 

The first step is to construct for two dimers a "bare" interaction potential, or vertex, A4, which is the sum of all 
irreducible diagrams, and then to build ladder diagrams from this vertex, in order to obtain an integral equation 
(see Fig. (S)). These irreducible diagrams are those ones which cannot be divided by a vertical line into two parts 
connected by two dimer lines. As it was pointed above the vertex A4 is given by the series shown on Fig^, since the 
diagrams on Fig^are by contrast reducible. Again we can eliminate Ts from our considerations and express A4 only 
in terms of T2. For this purpose we have to introduce a special vertex with two fermionic and one dimer incoming 
lines and two dimer outgoing lines Tai3{qi,q2',P2i P) (see Fig. I10|) . This vertex Tai3{qi,q2',P2, P) corresponds to the 
vertex A4 with one incoming dimer line being removed, in much the same way as ^{qi, q2]P2T P) and T/^{pi,p2 \ P) 
are related in Eq. ^ . The difference is that {qi ,q2;p2,P) is irreducible with respect to two dimer lines while 
$((7i, q2\P2, P) is not, just in the same way as T4^{pi,p2] P) and l^4^{pi,p2\ P) are related. The corresponding equation 
relating Taf3iqi,q2;P2,P) and A4{pi,p2;P) is: 

A4(pi,P2;P) = ^ xia,P)G{P+Pi~Q)G{Q)r^p{P + p,-Q,Q;p2,P). (Bl) 

Q; a, 13 

One can readily verify that the diagrammatic expansion for F shown on Fig. ^2 yields the same series as the 
one shown on Fig. for the vertex A4. The spin part of Ta,i3 has again the simple form Ta^/siqi, q2', P,P2) = 
x(a, /3)F(gi, q2',P2, P) and the function F(gi, q2]P2, P) obeys the following equation: 

r(gi, 92; P2, P) = -G{P -qi+ P2)G{P -q2- P2) - G{P -q2+ P2)G{P - qi - P2)- 

- G{Q)G{2P ~qi~q2-Q) [T2{2P - qi - Q)T{qi, Q;p2, P) + ^2(2^ - q2 - Q)F(Q, 92; P2, P)] ■ (B2) 
Q 
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FIG. 11: The graphic representation of the equation on the full vertex T{q\,q2;p2,P)- 



The sign minus in ljB2|) is a consequence of the anticommutativity of Fermi operators. It is clear that Eqs. ljBl|) and 
(|B2|I can be analytically integrated over the variable fi. Thus the s-wave component of the vertex r(gi, (72;?'2, -P) is 
a function of the absolute values of vectors |qi| and |q2|, the angle between them, the absolute value of vector |p2|, 
and the frequency u>2- The s-wave component of the sum of all irreducible diagrams A4(pi,p2; P) is a function of the 
absolute values of the vectors |pi| and |p2| and the frequencies wi and tJ2- 

The fully symmetrized vertex T4{pi,p2; P) of two-dimer scattering can be found from the solution of the following 
equation (see Fig. |^: 

T4(pi,P2; P) = A4(pi,P2; P) + IY. ^4(pi, g; P)T2(P + q)T2iP - q)T4q,p2;P). (B3) 

where A4(pi,p2; P) is the sum of all irreducible diagrams, P ± pi,2 — {—Eh ± wi,2, ipi.2} are 4-vectors of incoming 
(1) and outgoing (2) dimers in the center-of-mass system. 

Let us finally note that, equivalently to our above derivation, Eq. (|BllB3p can be also related to Eq. ||SJ| and Eq. 
© algebraically by simple formal operator manipulations. 
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